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Abstract. In this article we present an analytical cal- 
culation of the probability distribution of the magnifica- 
tion of distant sources due to weak gravitational lensing 
from non- linear scales. We use a realistic description of 
the non-linear density field, which has already been com- 
pared with numerical simulations of structure formation 
within hierarchical scenarios. Then, we can directly ex- 
press the probability distribution P{fi) of the magnifica- 
tion in terms of the probability distribution of the density 
contrast realized on non-linear scales (typical of galaxies) 
where the local slope of the initial linear power-spectrum 
is n = —2. We recover the behaviour seen by numerical 
simulations: P(/i) peaks at a value slightly smaller than 
the mean (/i) = 1 and it shows an extended large fi tail 
(as described in another article our predictions also show 
a good quantitative agreement with results from N-body 
simulations for a finite smoothing angle). 

Then, we study the effects of weak lensing on the 
derivation of the cosmological parameters from SNela. 
We show that the inaccuracy introduced by weak lens- 
ing is not negligible: Aflm ^ 0.3 for two observations at 
Zs — 0.5 and Zg = 1. However, observations can unam- 
biguously discriminate between Qm — 0.3 and flm = 1. 
Moreover, in the case of a low-density universe one can 
clearly distinguish an open model from a flat cosmology 
(besides, the error decreases as the number of observed 
SNela increases). Since distant sources are more likely to 
be "demagnified" the most probable value of the observed 
density parameter flm is slightly smaller than its actual 
value. On the other hand, one may obtain some valuable 
information on the properties of the underlying non-linear 
density field from the measure of weak lensing distortions. 



gravitational lensing 
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1. Introduction 



Using Type la supernovae as standard candles, it is pos- 
sible to derive the cosmological parameters and fl^ 
from the observed magnitude-redshift relation (e.g. Perl- 



mutter et al.l999). However, several sources of uncertainty 
(e.g. redshift evolution of the luminosity of SNela) can af- 
fect this method. In particular, the apparent magnitude 
of distant supernovae can be distorted by gravitational 
lensing by the density fluctuations along the line of sight. 
In view of the importance of the measure of Qm and rt\ 
from SNela, it is of interest to get a good estimate of 
the effect of weak gravitational lensing. Moreover, since 
the latter is directly linked to the matter distribution in 
the universe, one might use this effect to get some infor- 
mation on the large-scale structure of the universe itself. 
Note that contrary to the usual weak lensing statistics ob- 
tained when one considers the filtered distortions realized 
on large angular scales to probe the quasi-linear regime 
(e.g. Bernardeau et al.l997), the effects studied in this 
article which are relevant to SNela come from strongly 
non-linear scales 100 kpc). 

Several authors have already studied some aspects of 
weak gravitational lensing in this non-linear regime, by an- 
alytical means (e.g. Frieman 1997; Kantowski 1998; Met- 
calf 1999; Hui 1999) or numerical simulations (e.g. Wamb- 
sganss 1997; Jain et al.l999). Here we present an ana- 
lytical calculation of the probability distribution of the 
magnification due to weak lensing, from a model of the 
non-linear density field which has already been compared 
with numerical simulations of structure formation. Thus, 
we can directly express the properties of the weak lensing 
magnification in terms of the characteristics of the un- 
derlying density field. In particular, we recover the non- 
gaussian behaviour of the magnification as seen in previ- 
ous numerical studies. Then we show that the the fluctu- 
ations of the magnification {6fj, ~ 0.08 at Zg = 1) lead to 
a significant uncertainty for flm (Afim ^ 0.3 for two ob- 
servations at = 0.5 and z^ — I, but this error decreases 
as the number of observations increases). 

This article is organized as follows. In Sectj^ we recall 
the description we use for the density field. Next, in SectJ^ 
we derive the probability distribution of the magnifica- 
tion of distant sources by weak gravitational lensing. We 
present in Sect.^ the numerical results we obtain for three 
cosmologies (critical, open and low-density flat universes) 
as well as the dependence on the cosmological parameters 
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of the amplitude of the fluctuations of the magnification. 
Then, in Sect.| we compare our results with another ap- 
proach, used by several authors (e.g. Porciani & Madau 
1999), where the density field is described as a collection 
of smooth virialized halos. In particular, we point out the 
limitations of this method (it is restricted to large mag- 
nifications and it leads to some inconsistencies with some 
results from numerical simulations). Next, in Sect|6| we 
show how weak lensing can affect the measure of fim and 
r^A from SNela. Finally, in Sect .0 we present the bias due 
to weak lensing, that is the distortion of the luminosity 
function of SNela. 

2. Density contrast probability distribution 

In order to obtain the probability distribution of the mag- 
nification of distant sources by gravitational lensing we 
need the properties of the density field. Hence we briefiy 
recall here the formalism we use to characterize the density 
fluctuations. We shall use the same techniques to derive 
the statistics of the flux perturbations. It is convenient to 
express the probability distribution of the density contrast 
S = {p — p)/p at scale R and redshift z (here p{z) is the 
mean universe density) in terms of the many-body correla- 
tion functions ^p(ri, ...,rp). Thus we define the quantities 
(P>2): 



n — — ^ with = 



^2 



d^ri...d^i 



VP 



(1) 



where V = Att/3R^ is the volume of a spherical sphere of 
radius R. Next we introduce the generating function 



= 1 



(2) 



with Si — 1. Note that the parameters Sp can also be writ- 
ten in term of the cumulants {5^)c of the density contrast 
5 at scale R: 



(3) 



Then, one can show (White 1979; Bahan & Schaeffer 1989) 
that the probability distribution of the density contrast 5 
within spheres of size R is: 



P{5) 



dy 
2mS, 



(4) 



where we note ^2 The relation (y) provides the link 
of the density probability distribution with the correlation 
functions, hence with the cumulants of the density field. 
In the non-linear regime ^ ^ 1 it is convenient to define 
the function (inverse Laplace transform): 



hix) = - 



(5) 



which obeys: 

p>l : Sp^ xPh{x) dx , 81^82^1 (6) 

From very general considerations (Balian & Schaeffer 
1989) one expects the function (f{y) to behave as a power- 
law for large y: 

y +00 : ip{y) ^ a y^"" with < uj < 1 (7) 

and to display a singularity at a small negative value of y: 

y-^yt : f{y) = -a^ r{ujs) {y - y^)""' (8) 

where we neglected less singular terms. From this be- 
haviour of ip{y) we have: 



a; < 1 : h{x) ~ -^r—^ x'^"^ 
x:$>l : h{x) - a, e-^/^= 



(9) 



with Xs = l/\ys\- The interest of the function h{x) is that 
in the non-linear regime for large density contrasts the 
density probability distribution can be written as (Balian 
& Schaeffer 1989): 

? » 1 , (1 + <5) » e : P{S) = i h{x) (10) 



with: 



1 + 5 



(11) 



Thus, the density probability distribution P{5) shows a 

power-law behaviour from [\ + 5) ^ ^ up to 

(1 + (5) ~ XaS, with an exponential cutoff above Xs£,- The 
measure of P{5) in numerical simulations allows one to 
recover h{x) hence '^{y) since (|^) can be inverted as: 



(1 - e-'^y) h{x) dx 



(12) 



The fmiction h{x) has been measured in the non- 
linear regime for various power-spectra by several au- 
thors (Valageas et al.l999; Bouchet et al.l991; Colombi 
et al.l997; Munshi et al.l999). In particular, although 
Colombi et al.(1996) found a small scale-dependence other 
authors found that the numerical results were consis- 
tent with h{x) being scale-independent in the non-linear 
regime. Thus, in the following we shall use the scaling 
function h{x) obtained by Valageas et al.(1999) for any 
redshift. Note that it depends on the power-spectrum and 
it must be obtained from numerical simulations since there 
is no known method to derive analytically h{x). The scale- 
invariance of h{x), hence of the coefficients 5^, can be in- 
terpreted as evidence for the stable-clustering ansatz (Pee- 
bles 1980): 

ep(Ari, \vp- a) = a3(p-i) ^-7(p-i) ^^^^ ...^j-p) (13) 



p. Valageas: Weak gravitational lensing effects on tlie determination of f2,„ and Qa from SNela 



3 



which was studied in details in Balian & SchaefFer (1989). 
Here a{t) is the scale-factor and 7 is the (local) slope of 
the two-point correlation function. Note that for an initial 
linear power-spectrum which is a power-law P{k) oc we 
have if stable-clustering is valid: 



7 = 



3(3 + n) 



(14) 



The interest of the formulation (Q) is that once h{x), 
or ip{y), is known the density probability distribution 
P{S) can be obtained for any time and scale in the non- 
linear regime provided that one knows the behaviour of 
^ (and that indeed h{x) is scale-invariant). Note that a 
similar technique can be used in the quasi-linear regime 
with a different (p{y) obtained by perturbative calculations 
(Bernardeau 1994). 

Then, in order to obtain the properties of the non- 
linear density field we only need to model the evolution 
of the two-point correlation function ^2 , or of the power- 
spectrum P{k). To this order we use the fits given by 
Peacock & Dodds (1996) which give the non-linear power- 
spectrum P{k) from its linear counterpart P^ik). Note 
that this behaviour of P{k) is consistent with the stable- 
clustering ansatz. 

3. Magnification by weak lensing 

3.1. Definition 

As a photon travels from a distant source towards the ob- 
server its trajectory is deflected by density fluctuations 
along the light path. This produces an apparent displace- 
ment of the source as well as a distortion of the image. In 
particular, the convergence k (defined as the trace of the 
shear matrix) will magnify (or demagnify) the source as 
the cross section of the beam is decreased (or increased). 
One can show (Bernardeau et al.l997; Kaiser 1998) that 
the convergence is given by: 



3 /"^^ 
K = -flm J dx w{x,Xs) 5{x) 



with 



w{x,Xs 



Hi VjxMXs - x) 
-Dixs) 



(1 + ^) 



(15) 



(16) 



where x is the radial comoving coordinate (and Xs corre- 
sponds to the redshift Zs of the source): 



dx 



-n- dz 

no 



+ (1 - r!„, - r!A)(l + Z)2 + an(l + ZY 

while T) is defined by: 



sinh ( \/l — ilm — 



(17) 



(18) 



The relation (|l|) assumes that the metric perturbations 
(f) are small {(f) ^ 1) but the density fluctuations 5 can be 



large (Kaiser 1992). The magnification /i is linked to the 
convergence n and the intensity of the shear 7 by: 



1 



(19) 



Thus, for small values of n we have: 

K<1 : h=1 + 2k (20) 
and we write the flux perturbation 5^ — ^ ~ \ as: 

(21) 



5n = 3n,n dxw{x,Xs) S{x) 
'0 



Of course, because of flux conservation the mean shift of 
the flux over all lines of sight is zero: {5fi) = 0. We can see 
directly in (^ij) that there is a minimum value /^niin(^s) 
for the magnification of a source located at redshift Zg: 



lLmin{Zs) = 1 - 3r2,„i^s(Xs) 

with 



FsiXs 



dx w{x,Xs 



(22) 



(23) 



This corresponds to an "empty" beam between the source 
and the observer (S = —1 everywhere along the line of 
sight). 

3.2. Probability distribution P{^) 

Next we wish to obtain the probability distribution of the 
magnification from (pT|). To this order we simply need 
to derive the cumulants {Sfj,P)c. This will provide the pa- 
rameters Sfj,^p, similar to (^, and the generating function 
(Pf_i{y), similar to (^. However, it is convenient to intro- 
duce first a "reduced magnification" 77 and its variation Srj 
by: 



1 - 



, Sr] = rj — I 



(24) 



From ( |2l| ) and ( p4[ ) we obtain the cumulant of order p of 
6rj: 



dXi -E, KXi) 



(25) 



Since the correlation length (beyond which the many-body 
correlation functions are negligible) is much smaller than 
the Hubble scale c/H{z) (where H{z) is the Hubble con- 
stant at redshift z) we obtain: 



Jo 



F. 



00 P 



(26) 



.Xp; 
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From ( [T3I ) we see that the integral over the points Xi 
along the line of sight is dominated by the comoving scale 
Rc{z) such that the local slope of the two-point corre- 
lation function is —"/(Re) = —1. This corresponds to a 
local linear power-spectrum index n — —2, see (p^). In- 
deed, for realistic power-spectra like CDM the slope of 
^2{x) decreases from at small scales down to —4 at large 
scales. On the other hand, power-spectra which would be 
pure power-laws would lead to divergences. Since for the 
power-spectra we shall use Rc{z) corresponds to galac- 
tic scales (i?c(0) ^ 100 kpc) within the highly non- linear 
regime (^(i?c;0) ~ 2000), we must indeed use the non- 
linear many-body correlation functions Cp(0, X2j •■, XpS 2;). 
Thus we can use the scaling laws (|3|). Of course, at large 
redshifts z > 7 the scale Rc{z) will enter the linear regime, 
however since we shall restrict ourselves to smaller z we 
can always use ([T^). Now we must estimate the contribu- 
tion of the integral over the points Xi along the line of 
sight. Although the points Xi £^re drawn on a line and not 
within a sphere we shall use the approximation (compare 
with (0)): 

/oo P / poc \P— 1 

Yldxi ^p{o,x2,--,xp) - c^x6(o,x)) 

(27) 

where the coefficients Sp are obtained for n = —2 (i.e. 
7 = 1) from the numerical results described in Valageas 
et al.(1999). Thus, we define the quantity Ifj.{z) by: 



(28) 



/^(z) = / dx C2ix;z) 



It is convenient to express in terms of the non-linear 
power-spectrum which is directly provided by the fits ob- 
tained in Peacock & Dodds (1996). Thus, using the Fourier 
transform: 



6(x)= / d^fce 



P^k)= r^A2(fc)£m(M (29) 
Jo ^ "--^ 
where we defined the power-spectrum P{k) by: 

((5(ki)<5(k2)) = P(fci) SdO^i + k2) , A2(A:) = 47rfc3p(fc) 

we obtain (note that 2; and k are comoving coordinates): 

dk A2(fc) 



(30) 



Next, using the approximation ( |27| ) in the relation (g6|) we 
get: 



Sp If,{z) 



p-i 



(31) 



In particular, this yields for the variances — {^tf') and 
Cm = {^^^^) of the fluctuations of the amplification of a 
source located at redshift Zc: 



I,.. 



= (Sn^f J dx (32) 



Note that our result ( p^ ) for the rms fluctuation does 
not use the approximation (|2^). Moreover, the expressions 
( ^2[ ) clearly show that most of the contribution to weak 
gravitational lensing effects comes from the scale Rc{z) 
where A'^(k)/k is maximum. This corresponds to 7 = 1 
or n' = —2 where n' is the slope of the non-linear power- 
spectrum. From ( p^ ) we see that it also corresponds to 
n = —2 where n is the slope of the linear power-spectrum. 
Then, from (|3^) we can define the generating function 
(Prfiy), similar to (|), by: 



This yields: 



(-l)f-i {57^P), 



p=2 



fP-1 



"Xs ^ / w I 



(33) 



(34) 



where f is the generating function of the density contrast 
defined in (^). Next, in a fashion similar to (^) since {Srf) = 
0, we obtain for the probability distribution of 77: 



P(ry) 



dy 



:,[VV-Vr,{y)]/ir, 



and the p.d.f. P{fi) is given by: 
1 



1 A^min 

The expansion of ifriiv) in y = verifies: 
1/ 

f^iy) = y-'— + ... 

which implies, as it must, that: 

(^) = 1 and {{t^ ^ {r^)f) ^ 



(35) 



(36) 



(37) 



(38) 



The relations (|34| ) and ( p5|) provide the probability distri- 
bution of the magnification of a source located at redshift 
Zs- Note that the generating function fniy) depends on 
Zg. We can also check in (35), using (|3^), that P(?7) = 
for 77 < as it should (since for 77 < we can push the 
integration path in ( p5| ) towards the right, Re(y) -l-oo, 
where the exponential vanishes). Thus, the approximation 
has preserved the fact that P{fJ,) — for fi < /imin- 
In a fashion similar to (^ we can define the function: 



hrj(x) 



-\-ioo 



dy 

2TTi 



dx¥ 



1 



(39) 



with: 

Ts{x) 



(40) 
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The function hri{x) satisfies: 



with 

p>2 



{61]' 



X^h,-i{x) dx , Sri^i — Sri,2 = 1 (41) 



(42) 



p-1 



From ( |39| ) and (^ we can obtain the parameters which 
govern the asymptotic behaviour of hrj{x): 



ri.s — s T^sijCc) 
1/2 



(43) 



t:{xc) 



where < Xc < Xs is the point where T's(x) is maximum. 
Thus, we see that the slope of /i,,(a;) at low x is the same 
as h{x) while the exponential cutoff is slightly modified 
{Tsixc) ^ !)■ From (||) and (|3|) we see that T^ix) ^ 1 
so we are led to the approximation: 



hrj{x) = h{x) and (prj{y) = f{y) 

This leads to the approximate p.d.f. for 77: 



-\-ioo 



dy 



(44) 



(45) 



which is still properly normahzed and satisfies P{t]) — 
for 1] < 0. The practical advantage of the approximation 
( ^ ) is that to compute the p.d.f. P{r]) we can directly 
use the functions h{x) and ip{y) obtained in numerical 
simulations. Note that formally P{r]) can be expressed in 
terms of the probability distribution P{5) of the density 
contrast 5 at the scale such that n = — 2 as: 



(46) 



using (||) . Note however that even for < 1 the probabil- 
ity distribution P{S) used in ( p6| ) corresponds to the highly 
non-linear regime. Thus in the case < 1 the probability 
distribution P{6) in ( ^ ) is not the one measured at the 
time such that ^ = Cij- Nevertheless, the expressions ( p5| ) 
and ( p6| ) clearly show that the measure of the p.d.f. P{ij), 
or of P(k), provides a direct estimate of the p.d.f. P[5) of 
the density contrast through Lp{y) (see also the study of 
P{k) in Valageas 1999b). 

As can be seen from ( |3^ or (^) the probability dis- 
tributions P{rf) and P{n) are non-gaussian. Indeed, they 
are strictly zero for r/ < and fi < //min, they show an ex- 
ponential cutoff ^ Q-'n/{xr,.sir,) (■qj. ^ Q-^J■/{x■n,s{^~^J.miT,)il^)^ 

for large 77 (or large ^), and they have their maximum at 
a value smaller than the mean {rf) = (fi) = 1. 



We compare in details the predictions of our approach 
with available numerical results from N-body simulations 
(Jain et al.l999) for the convergence k smoothed on small 
angular scales {6 ~ 1') in Valageas (1999b). This com- 
parison shows that our approach (which can be extended 
to finite smoothing windows in a straightforward fashion) 
provides very good results for all three cosmologies we 
consider here (e.g., see Fig. 4 and Fig. 5 in Valageas 1999b). 
Moreover, we show in that paper that the approximation 
(j^) gives reasonable results which are quite close to the 
more accurate expression (jsj). Indeed, we shall check be- 
low that the correction to the third moment (for instance) 
due to (|4|) is quite small and well within the errorbars of 
the value of S3 obtained by counts-in-cells calculations 
from numerical simulations. Note that the evaluation of 
the probability distribution of the magnification fi by ray 
tracing through N-body simulations would of course suffer 
from the same uncertainty, which affects the moments of 
the density distribution itself (because of numerical inac- 
curacy). 

3.3. Scaling functions ifiy) and h{x) 

In order to perform numerical calculations, we need to 
choose the function h{x), or equivalently the parameters 
Sp defined in (0). As we explained above, most of the 
contributions to the weak lensing effects come from the 
scales where the local slope of the linear power-spectrum 
is n = — 2, see also Fig. 2 in Valageas (1999b). Hence we use 
in the following the scaling function h{x) obtained from 
numerical simulations by Valageas et al.(1999) for the case 
of a critical universe with an initial linear power-spectrum 
which is a power-law n = — 2: 



h{x) 



with: 



a{l-Lu) x"^ 



T{uj) {l + hxY 



exp(-a;/xs) 



(47) 



a = 1.71 , u = 0.i 



13 , 6 = 5 and c = 0.6 



hence ujs = uj — 1 — c = —1.3 

In fact, the curvature of the CDM power-spectrum may 
slightly change the parameters Sp from the value they 
would have for a pure power-law P{k). However, in order 
to improve meaningfully this approximation one would 
need to measure the parameters Sp realized on a line 
rather than in a sphere, see (p7|). Thus, we think our ap- 
proach is the best analytical tool one can currently build. 
The scaling function h{x) shown in ( ^7| ) defines the gen- 
erating function (pijj) through (^2|). In particular, one ob- 
tains (see Gradshteyn & Ryzhik 1965, §9.211, p.l058): 



ip'{y) = a(l - uj)b^'^ ip[uj,2 + uj. 



y-ys 



(48) 



6 



P. Valageas: Weak gravitational lensing effects on tlie determination of f2,„ and Qa from SNela 



and 



Table 1. Cosmological models 



(p{y) = a 6^ 



4' a; — 1 , 1 + 



y-Vs 



(49) 



i/; { Lli - 1,1 + ujs 



b 



where tp is Kummer's function which can be expressed in 
terms of the difference between two confluent hypergeo- 
metric functions iFi. Next, to obtain P{ri) from (45) one 
simply needs to perform an integration in the complex 
plane. In order to make the integral ( ^ ) converge suffi- 
ciently fast, it is convenient to define the integration path 
by the constraint Tm[rjy — ^p{y)\ = 0. However, in prac- 
tice it is sufficient to use lm[r]y — a{y — ys)^""^] = (we 
replace f{y) by a power-law with the right behaviour for 
y — *■ +00 and the right location of the singularity at y^) 
which gives for the integration path: 



with 



a sin(l — uj)6 
T] sin 9 



(50) 



Note that it is better to define (p{y) from h{x) rather than 
trying to use a fit for (p{y) itself. Indeed, from (|l]) and (||) 
we see that: 



p>l 



-1)P- 



-V^''^(o) 



Sp>0 



(51) 



and moreover, using Schwarz' inequality and the scalar 
product {f\g) — J f{x)g{x)xh{x)dx, one can see that the 



coefficients Sp must obey: 



Sp-q 



> Si 



(52) 



These constraints are automatically verified if one defines 
ip{y) from h{x). If one uses a fit for ^(y) which does not 
obey these constraints one may get negative probabilities 
(since in this case h{x) has to be negative in some range). 

The fit (^) for the function h{x) was obtained for a 
critical universe. In the case of a low density universe, we 
use the same function although there are no numerical re- 
sults available (from counts-in-cells statistics) to validate 
(or invalidate) this choice. However, we note that the fact 
the dependence on cosmology of the two-point correlation 
function is accurately given by the simple term described 
in Peacock & Dodds (1996) suggests that the structure of 
the non-linear clustering pattern is the same for low Q„i 
as for a critical universe, once the effect of the slow-down 
of the linear growth factor is taken into account. Indeed, 
most highly non-linear structures formed when the uni- 
verse was close to critical {flm{z) ~ 1) since at later times 
the slow-down of the linear growth factor prevents addi- 
tional new structures to form. Of course, this break in the 
hierarchy of scales which successively turn non-linear may 
also lead with time to some difference with the case of 
a critical universe (at least for the scales which were the 
last to collapse) . Detailed numerical studies are needed to 
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investigate more precisely this point. However, for reason- 
able cosmologies ilm ^0.1 our model provides a good ap- 
proximation as shown by a direct comparison of P{k) with 
results from N-body simulations, as described in Valageas 
(1999b). 

Note that, except for this possible dependence of the 
parameters Sp, all our results are valid for any realistic 
power-spectrum such that n < —2 on small scales, n > —2 
on large scales and the scale where n = — 2 is non-linear. In 
particular, note that for such power-spectra all moments 
of the convergence k and of the magnification ^ converge: 
there is no need to introduce a cutoff at small scales. 

4. Dependence on cosmology and redshift 

Using the results obtained in the previous sections we can 
compute the probability distribution P(/i) of the magni- 
fication for various cosmologies. We shall mainly consider 
three cases: a CDM critical density universe (SCDM), a 
low-density open universe (OCDM) and a low-density fiat 
universe with a non-zero cosmological constant (ACDM). 
The main cosmological parameters of these scenarios are 
described in Tab.|l|. Here T is the usual shape parameter of 
the power-spectrum and we use the fit given by Bardeen 
et al.(1986) for P{k). These cosmologies are those we use 
in Valageas (1999b) to compare our predictions with avail- 
able results from N-body simulations. The reader is refered 
to that paper for a detailed discussion of the accuracy of 
our approach and for an extension of our model to finite 
smoothing windows. 

4-1. Fluctuations of the magnification 

First, we present in Fig.^the redshift evolution of the am- 
plitude of the fiuctuations of the magnification /x of distant 
sources located at Zg . The increase with Zg of the interval 
(1 — Umin) between the mean {/i) = 1 and the minimum 
value /imin is due to the more extended line of sight which 
gives more room for the weak lensing effects. This devia- 
tion is smaller for the low-density universes than for the 
critical case because of the factor which enters (p2[): 
the difference of matter between the mean and is smaller 
as it is proportional to f2m at low z. It is larger for the 
flat cosmology (with the same J7m) because of the detailed 
dependence on 17 a of the factor Fg (see also Bernardeau et 
al.l997): indeed the distances x and V are larger at fixed z 
which gives more room for weak lensing effects at a given 
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Fig. 1. The fluctuations of the magnification /i for critical, 
open and low-density flat universes, for a source located at 
redshift Zg . The solid lines show 1 — fimin where jimin given 



by (22) is the minimum value of the magnification. The 
dashed lines present the variance y/£,fj, of the magnification, 
from (|32[ ). The dotted lines show the variance y/^rj of 77 
from (|32|). 

Zg. On the other hand, the variance -^/^ of the "reduced 
magnification" rj decreases at larger redshift because the 
integration in over the successive "slabs" of matter 
leads to be "closer" to a gaussian, in a fashion similar 
to the central limit theorem (although the latter does not 
apply here since the probability distribution of the magni- 
fication due to each slab evolves with z and Zg). It diverges 
for Zg where the number of (highly non-linear) den- 
sity fluctuations which intersect the line of sight declines. 



Note however that for large fi the probability distribution 
P(/x) is always very different from a gaussian, whatever 
the value of and since it shows a simple exponential 
cutoff rather than a gaussian cutoff. As a consequence, at 
low Zg the variance of the magnification /i becomes of 
the order of, and even larger than, (1 — ^min) (thus P(/i) 
is strongly non-gaussian and sharply peaked close to the 
minimum fj-min) while at large redshift becomes sig- 
nificantly smaller than (1 — ^min) (thus the peak of P{fJ.) 
gets closer to the mean (/i) = 1). In particular, we can 
check from (||) and (H) that: 



: (1 



Mil 



1) oc , £,fj, oc Zg and cx Zg 



and: 



00 : is finite and cx (1 + Zg) 



(53) 



(54) 



We show in Valageas (1999b) that our results agree with 
the values obtained by Jain et al.(1999) using ray tracing 
through N-body simulations, for smoothing angles 9 ^ 
0.1'. In fact, since our prediction for the variance only 
relies on the weak lensing approximation and on the fits 
for the non-linear power-spectrum given by Peacock & 
Dodds (1996) this agreement mainly shows that both sets 
of simulations are consistent. As seen in Fig.l in Valageas 
(1999b), the variance obtained without smoothing {9 = 
0), which we consider here, is slightly larger than the value 
reached at 6* = 0.1' (since a finite smoothing removes the 
power from small scales). For instance, for the SCDM case 
we get ~ 0.1 for = (as in Fig.|l| here) and ~ 0.08 
for 9 = 0.1'. 

4-2. Probability distributions 

We display in Fig.|^ the probability distribution P{r]) of 
the "reduced magnification" 77, from (^). As explained in 
Sect.^ it is strongly non-gaussian with a maximum below 
the mean (77) = 1 and a power-law tail followed by a simple 
exponential cutoff at large 77. Its is more strongly peaked 
around its maximum which is closer to the mean (77) = 1 
for the ACDM model, and even more for the SCDM sce- 
nario, following the decrease of shown in Fig.|l| (but 
of course this depends on the cosmological parameters one 



chooses). Note that even for very small variance -^^^ of 
the magnification 0.1, see Figjl]) the probability 

distribution of the magnification displays clear deviations 
from gaussianity, as explained above. In particular, the 
normalized magnification 77 clearly shows the shape of the 
probability distribution, which appears squeezed towards 
the mean (/i) — 1 when displayed as a function of /x. 

This is shown in Fig.|^ where one can check that at 
low redshift P{n) tends to a Dirac 5d(/x — 1). Due to the 
smaller value of (1 — fimin) and for low density uni- 
verses P{p) is much more sharply peaked around its max- 
imum than for a critical cosmology. One can clearly see 
the asymmetry due to the lower cutoff at fj-min and the 
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Fig. 2. The probability distribution P{r]) of the "reduced 
magnification" 77, from (p5|). The solid lines correspond to 
Zs — I and the dashed lines to Zg — 0.5. 

extended large fi tail, as well as the shift of the maxi- 
mum of P{fJ.) below the mean (/i) = 1. This is similar to 
the behaviour obtained from numerical simulations (e.g. 
Wambsganss et al.l997). We display below in Fig.|^ the 
curve log[P(/i)] which shows even more clearly the non- 
gaussian features of -P(/i). For the critical density uni- 
verse we also show in Fig.|| the gaussians (dotted lines) 
which correspond to the same variance We can see 
that the probability distribution P{p) is indeed very dif- 
ferent from a gaussian, at both redshifts. Thus, it would 
be quite meaningless to model P{fJ.) as a gaussian. For the 
sake of completeness, we note here how one can obtain the 



Fig. 3. The probability distribution P{fJ,) of the magnifi- 
cation ^, from ( |45| ) and (|36|). The solid lines correspond 
to Zs = 1 and the dashed lines to Zg — 0.5. For the crit- 
ical density universe we also show the gaussians (dotted 
lines) which correspond to the same variance ^^j. These 
gaussians have a peak at the mean (/i) = 1 and the larger 
redshift Zs — 1 corresponds to the larger variance and to 
the lower height of the maximum. 

the peak of the probability distributions P{ri) and P{n)- 
From (^) and (^ one can show that for large variance 
S,n the location Tymax of the maximum of P{'q) is given by 
(Bahan & Schaeffer 1989; Valageas & Schaeffer 1997): 

» 1 : 77max ^ ^--/(i-) (55) 

which implies: 



approximate behaviour of the locations 77max and /Xmax of > 1 : ^i^^^ ~ -j- (1 - 



S.J) 



(56) 
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Thus, for large the peak of the probability distribution 
of the magnification gets very close to the minima rj — 
and /Xmin- For smaller values of the variance one may 
write an Edgeworth expansion of the probability distribu- 
tion P(r]) (e.g. Bernardeau & Kofman 1995): 



1 



1 



S3 

' 6 



(57) 



where He^^i') = — "iv is the Hermite polynomial of 
order 3, v — (ry — l)/^;, and Un — This gives for the 

location of the peak of the probability distributions: 



< 1 : ?7n 



and 



< 1 : ^i, 



l-^(l-/i„ 



(58) 



(59) 



Thus, for both large and small redshift, that is for small 
and large <^^, the peak /imax tends to the mean 1 (but 
for different reasons). Hence the deviation (1 — /imax) is 
maximum for an intermediate redshift of order unity. 

From the probability distribution P(/i) one can obtain 
intervals of confidence for the magnification /i. Thus, for 
< < 1 we define , as the minimal interval 
(i.e. with the smallest length) such that /x € [/i_, with 
probability V . These intervals contain the location /imax 
of the peak of P(/t) and obey: 



P(/i_) = P(/^^ 



P(/i_ < /i < /i+) 



(60) 



They are displayed in Figji] for the redshifts Zg = 0.5 and 
Zs = 1: any horizontal line of ordinate V intersects the 
curves at the points /i_ and /i+. For V the length 
of the interval goes to as /i_ and fi-^- tend to /imax- For 
■p = 1 we have /t_ = /imin and /t+ = cxd. Thus, Fig.^ 
clearly shows the range of fi one can expect, as well as the 
asymmetry of the underlying probability distribution. In 
particular, the large ^ tail of P(/t) is clearly seen. More- 
over, one can check that there is a non-negligible proba- 
bility to have /i < 1 on a line-of-sight (/i+ < !)• 

4-3. Skewness 

Finally, in Fig.^ we display the third moment of the prob- 
ability distribution of the magnification. From (0) the pa- 
rameter S3 (skewness) is given by: 



r 



(61) 



since {5^) = {^r)c- It measures the third moment of 
the probability distribution of the density contrast (5i? 
realized in spherical cells of radius R. As explained in 
Sect.^ it is constant with time in the non- linear regime 
at the scale Rc{z) such that the local slope of the linear 
power-spectrum is n = —2. From the results obtained by 



0.6 



a, 

V 

=L 0-4 





Fig. 4. Intervals of confidence for the magnification fi. The 
curves show the probability P{fJ.- < 1^ < 11+) that the 
magnification fj, is within the minimal interval [/i_,/i+]. 
The solid lines correspond to = 1 and the dashed lines 
to Zs = 0.5. 

Valageas et al.( 1999) from numerical simulations we have 
53 ~ 10.7. From (45), 53 also measures the third order 
moment of r/: 



t2 - 

From the change of variable ( p^ we also obtain: 



s 



(62) 



(63) 
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Fig. 5. The parameters which define the third moment 
of the probabihty distribution for the magnification fi 
(5*^,3 and £^,3) and the "reduced magnification" r/ or 
the density contrast {S3). The sohd fines show the value 
of log(S'^,3) and log(S^,3) obtained from ( |6^ ) while the 
dashed curves correspond to the approximation (|63|). At 
larger redshift P{n) gets closer to a gaussian so that S'^,3 
and S^,3 decrease while S3 is independent of z^. 

This clearly shows the dependence on cosmology of the 
coefficients S',i,p. Moreover, from ( |3l| ) we can also derive 
the parameters without the approximation (^J). This 
leads to: 



Sfl,p (3^^771 -^s) ^ Sp 



dx 



p-i 



Finally, we also define the parameter S^.3 by: 



^M:3 — 



1/2 



(65) 



(64) 



which may be seen as a more convenient measure of the 
deviation of P(/i) from a gaussian. We can check in Fig|| 
that the error introduced by the approximation (Q) is 
quite small. In particular, it is negligible as compared to 
the inaccuracy due to the results of numerical simulations 
of structure formation. Indeed, while Valageas et al.(1999) 
get 53 ~ 10.7, Colombi et al.(1997) obtain S3 ~ 10.2 and 
Munshi et al.(1999) find 5*3 ~ 6. Thus, the approximation 
(Q is quite sufficient in view of the accuracy of the scaling 
function h(x) obtained from numerical results. However, 
one should note that the functions h(x) measured in simu- 
lations are reasonably close in the range where they have 
been tested against numerical data. Indeed, the uncer- 
tainty which affects the parameters Sp (and increasingly 
so for large p) comes from the large density tail of the 
probability distribution of the density contrast while the 
behaviour of x^h{x) around its maximum (i.e. at a; '--^ 1) 
is fairly well constrained, see Valageas et al.(1999) for a 
detailed discussion. Note that the measure of S^^p from 
ray-tracing in N-body simulations would of course bear 
the same inaccuracy, which is due to the dispersion of the 
properties of the non-linear density field itself obtained 
from different numerical simulations. 

However, we can note that at = 1 we get S'^_3 = 49, 
131 and 95 for the critical density, open and low-density 
flat universes, while Jain et al.(1999) obtain 5^,3 — 43, 
107 and 75 (we have 5^,3 = 8^,3/2). For the critical 
density universe both values agree quite well, while for 
the low-density universes our values are somewhat larger 
than those obtained by these authors (although they show 
the same trends). This could be due to the slow rise of 
the skewness with smaller smoothing angle (due to the 
fact that S3 is larger for non-linear scales than for quasi- 
linear scales, see Colombi et al.l997). Indeed, as seen in 
Jain et al.(1999) the skewness may not have reached its 
asymptotic value at 6* = 1' yet. Moreover, as discussed 
in Valageas (1999b) the errorbars on the measure of 5^_3 
from numerical simulations may be larger than the disper- 
sion of the estimator used to compute the skewness may 
suggests because two p.d.f. with a rather different skew- 
ness can still agree very well, as shown by the good agree- 
ment of our prediction for P{k) with the results from these 
N-body simulations. In particular, although this compar- 
ison may suggest a small dependence of S3 on flm the 
dispersion of numerical results (which provide values for 
53 which can vary by a factor 1.8 as seen above) prevents 
us from drawing definite conclusions. Of course, it would 
be interesting to measure both S3 and S'^,3 in the same 
numerical simulation to directly check the accuracy of our 
relation (|64|). On the other hand, the main improvement 
to our calculation would be to directly measure (or obtain 
from first principles) the parameters Sp defined by the in- 
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tegrals on the line-of-sight of the many-body correlation 
functions, i.e. the l.h.s. terms in (p7[ ), rather than the ra- 
tios of the averages defined in (|l|). Then one could still 
apply our method and simply use the new function ip{y) 
(or h{x)) defined by these new parameters Sp as in (^). 

As explained above, we can check in Fig.^ that at 
larger redshift the probability distribution of the magnifi- 
cation becomes "closer" to a gaussian as the parameters 
^1-1,3 and decrease. However, it is interesting to note 
that these parameters can be fairly large (higher than 6*3) 
and even diverge for z ^ 0. Thus, at low z the probability 
distribution of the magnification is strongly non-gaussian. 
In particular, we obtain: 

Zs ^ : Sfj_^3 cx z^^ and S^_3 oc zj^/^ (66) 
and 

Zs 00 : S'p^a oc (1 + Zs)"^ and S^^a oc (1 + z^)"^ (67) 

Thus, if the intrinsic magnitude dispersion of the sources 
is sufficiently small, one might observe these non-gaussian 
features and check that they agree with the usual mod- 
els of the density field. Moreover, from (^) we see that 
one could get an estimate of the properties of the under- 
lying density field itself (e.g. its skewness 53) from P(/i). 
However, as seen in Fig.yjthe width of the probability dis- 
tribution P(/x) is quite small at low z which would make 
such a study rather difficult, so that intermediate redshifts 
Zs ~ 1 may provide better results. 

4-.4- Influence of cosmological parameters 

We show in Fig.|| the dependence on Qm and fl\ of the 
fluctuations of the magnification n oi a source located at 
redshift Zg = 0.5 and Zs = 1. We vary the normalization 

— 1/2 

as of the power-spectrum with flm as as oc ilm which 
roughly accounts for the change of as needed in order to 
reproduce as well as possible large-scale structure observa- 
tions (abundance of clusters, velocity fields) with different 
cosmologies (of course, some cosmologies match such ob- 
servations better than others, whatever their choice of cts). 
Note that the shape of P{k) also depends on il.,„ through 
r = rimh. We can see in Fig.^ that (1 — ^min) increases 
for larger Qm- This is due to the factor Qm which appears 
in (p^). It translates the fact that for higher flm there is 
more matter in the universe hence there is more room for 
deviations from the mean (/i) = 1, for instance in the case 
where 5 = — 1 everywhere along the line of sight ("empty 
beam") which leads to fimm- In particular, we have: 

n„,^o : (1 

Mmin) OC ^jyi (^^) 

for both fiat and open universes, since for an empty uni- 
verse all beams are empty. The variances and 
increase for ilm because of the rise of the two-point 
correlation function which compensates the slower growth 
of the linear growth factor, see Peacock & Dodds (1996). 
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Fig. 6. The dependence on flm and fl^ of the fluctua- 
tions of the magnification /i of a source located at redshift 
Zs — 0.5 and Zg = 1. The solid lines show log(l — /imin)j 
the dashed lines present the variance log(-y<^^) of the mag- 
nification and the dotted lines show the variance log(-y^,,) 
of r]. The higher redshift corresponds to smaller and 
larger and (1 - Umin)- 

The amplitude of the fiuctuations of the density contrast, 
hence of the magnification fi, is smaller for a fiat uni- 
verse than for the open case with the same flm because 
of the detailed form of this linear growth factor. Overall, 
the variation of the quantities (1 — /imin), \/$n and 
is rather small over the usual range of cosmological pa- 
rameters 0.2 < n„i < 1. Note that the moments of the 
probability distribution P{ii) also depend on the cosmo- 
logical parameters and JIa (mainly on O^), as can be 
seen from (|63|). In particular, we have: 

5^,3 ^ — ^ (69) 

SO that the variation of 5^1,3 can be directly seen in Fig|| 
from the evolution of (1 — /imin)- 

5. Galactic halos 

In the previous sections, we have obtained the probability 
distribution of the magnification due to weak lensing by 
assuming the density field follows a specific scaling model 
described by (|l^), which is consistent with numerical re- 
sults (Valageas et al.I999; Munshi et al.l999; Colombi et 
al.l997). More precisely, although the parameters Sp and 
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the function h{x) may show some slow scale-dependence 
through the variation of the local power-spectrum index n 
they should not explicitely depend on scale (i.e. on cr^) so 
that the approximation ( p7| ) is valid. In turn, such a pic- 
ture of the non-linear density field implies that non-linear 
mass condensations can be described as an infinite hier- 
archy of increasingly small and overdense substructures 
(which follow the rise at small scales of as discussed in 
Valageas et al.(1999) and Valageas (1999a). On the other 
hand, it is customary to model collapsed objects as smooth 
halos with a density profile which follows a power-law or 
a curved profile as those obtained by Navarro et al.(1996) 
or Hernquist (1990). Moreover, within galactic cores where 
baryonic matter plays a dominant role one may expect a 
smooth density profile. Hence in this section we consider 
the case where collapsed halos are described by a power- 
law density profile: 



p(r) 



(1 



(70) 



which defines the slope /3 (with < /3 < 3). The factor 
(1 -I- Ac) is the overdensity at the virial radius R of the 
object. Recently a similar description was used by Porciani 
& Madau (1999), with f3 = 2 (isothermal sphere), with the 
Press- Schechter mass function (Press & Schechter 1974) 
of just-virialized objects. Here we intend to compare the 
results obtained from such a model to those we obtained 
in the previous sections. 

First, we note that such a description of the density 
field leads to a power-law tail at high /i for the proba- 
bility distribution of the magnification. Indeed, as shown 
in Valageas (1999a) the p— point correlation functions im- 
plied by (|7^) follow in the highly non-linear regime the 
power-law behaviour: 

p>2 : Cp(^) « R^'^" with 7p = p/3 - 3 (71) 
From (|26| ) we sec that: 

{Sfi ^ (Srj ^ RP-^ ^p - (72) 

where R is the upper or lower cutoff of the integrals in (p6|). 
Thus, if /3 > 1 and there is no lower cutoff (which would 
correspond to the radius of galactic cores where (3 < 1) 
the cumulants {6fj,P)c and the moments (/i^) diverge for 
p > 2/{j3 — 1). This is inconsistent with an exponential 
cutoff for P{ii) while it is the natural outcome of a power- 
law tail given by: 



large /i : P(/i) ^ /i 



-l-2/(,3-l) 



(73) 



We also note from the definition (|l|) that the description 
( [TOI ) means that the parameters Sp are scale-dependent: 



Sp oc i?-(P-2)(3-/3) 



(74) 



For reasonable values of /3 (/3 < 2) this implies a strong 
growth of the coefficients Sp [p > 3) at small scales which 
contradicts numerical results as discussed in Valageas 



(1999a). Note on the other hand that (3 < 2 leads to a 
small slope 7 < 1 for the two-point correlation function. 
Thus, we think this is not a good model for the dark mat- 
ter density field. However, we briefly present below the 
results obtained by such a description as it may be valid 
for baryonic cores of collapsed objects and it allows us to 
compare our previous results with other approaches. We 
restrict ourselves to 1 < /? < 3. 

Following the method used in Sect.^, in order to de- 
rive the probability distribution P{p) we first obtain the 
cumulants {6fi p)c- We divide the line-of-sight from the ob- 
server to the redshift Zg of the source into small physical 
elements Ak so that they contain at most one halo and 
from (|2^) we write for each realization the flux perturba- 
tion as: 



(75) 



where fii^a = 1 (resp. hi_a = 0) if there is (resp. there 
is not) a halo of type a within the length element i while 
Sfii^a is the contribution to the magnification of the source 
by a halo a at redshift Zi. From (Ell) we write: 



(76) 



S^li,a = S^m (1 + Z)'W{X, Xs) -f/3(l + A^) b ( — 



where: 



In ^2 



du 



(l+u2)/3/2 



(77) 



comes from the integration along the linc-of-sight through 
the halo. Here we made the approximation that the impact 
parameter b of the line-of-sight is much smaller than the 
radius R of the halo. The factor {1 + z) comes from the 
fact that here b and R (and I below) are physical scales 
while X is a comoving scale. Note that the index a denotes 
both the mass (or radius) of the halos and the impact 
parameter. Then, if we neglect the correlations between 
the collapsed objects we obtain in the continuous limit, 
using fii^a ^ 

{6f, P)c^ dl — db 27r6 (1 + z)3 ^^(M) 



M 



SQrn (1 + Z)w{x, Xs) Ipil + Ac) b 



R 



-13' 



(78) 



where we defined the comoving mass function of halos of 
mass M to M + dM as r]{M)dM/M. We also used: 



Ah db 2TTb (1 + z)^ r](M)^ 

M 



(79) 



Some of the integrals over the impact parameter diverge 
for P > 1 but we could add for these intermediate steps 
of the calculation an ad-hoc cutoff which we would put to 
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later on. In a fashion similar to (^) we can define the 
generating function ^pp^^i,{y) and we get: 



yA_ 



(80) 



with: 

A = (1 + z)w{x, Xs) + A,) b (^^ 

B = {1 + zf 2-Kh r]{M) 
and the variance of the magnification is: 



(81) 



= / dx 
Jo Jo 

Next, we define: 
(rih) ^ dx 



dM 



dM 
IvT 



dhB A" 



dbB 



(82) 



1 



and: 



dx 



dM 



(83) 



dbB A 



M/3 — A* ^ Ai/3,min 



(84) 



Thus, (n/i) is the mean number of halos along the line of 
sight while /U/j^min is the minimum value of the magnifica- 
tion Indeed, using (^ we obtain: 



P{pip) = e 



+ZOO 



^ W+fdX J^f dbBe- 



yA 



2ni 



(85) 



which shows that P{^i3) = for jip < 0. For large fj, such 
that P(///3)e<"''> < 1 we can develop the exponential and 
only keep the first two terms which leads for up > to: 

Pi^ip) ^ e-^"") JdxI^JdbBSnif^-A) (86) 

where i5d is Dirac's function. Of course, we could have 
written (^) directly, as it corresponds to the case where 
there is only one object along the line of sight. In practice 
{rih) ^ 1 and we obtain after integration over b: 



2tt 



dxl^R'viM) 



(87) 



x(l + z)2[3r!„, {l + z)w{x,Xs) W + ^c) R] 



We can check that we recover the power-law behaviour of 
([73|). Next, in order to compare with (45) we need to es- 
timate the halo mass function rj{M)dM / M . As discussed 
in details in Valageas & Schaeffer (1997) and Valageas & 
Schaeffer (1999), from the description of the non-linear 



density field presented in Sect.^ one can write the comov- 
ing mass function of halos defined by the density contrast 
A(A'f, z) (which may depend both on M and z) as: 



/ , ^\ dM Or. r, dx 
' M M ^ ' X 



1 + A(Af, z) 



where pg is the present mean universe density and h{x) 
is the function obtained from counts-in-cells statistics de- 
fined in d). A comparison of (|8|) with numerical results is 
described in Valageas et al.(1999). In order to get a direct 
comparison with (Ba) we now make the approximation: 



i{R;z) 



(1 + ^)" 



(89) 



This is possible because the halos we consider here cor- 
respond to galactic masses hence to the scale Rc{z) as 
noticed in Sect.^, where the local index of the initial lin- 
ear power-spectrum is n ~ —2. Here Rco = -Rc(O) and 
Cco — i{RcO', 0). Thus we eventually get: 



3(3»^J/3)~ 
2(/3-l) 



S2 {£,coRco) 



3-/9 



(90) 



20-3) 

dx {l + z)^^ w{x,Xs)' 



where we defined: 

/>oo 

= J xi^ h{x) dx (91) 
as in (^. In a similar fashion, we obtain from (|8^): 



3« r, p ^4(/3-i)fcr„ 



4r(cj) 



(92) 



x('f^) /dx(l + A.)^(l-fz)^ 



where the parameters a and oj describe the small x be- 
haviour of h{x), see (^, Ac ^ 177 is the density contrast 
of the smallest collapsed halos we consider at the thresh- 
old Tmin and /i™ is the mean molecular weight. Indeed, 
the number of halos along the line of sight is dominated by 
the contribution of the smallest objects (the multiplicity 
function diverges at small mass if there is no cutoff). We 
use Train — 3 10^ K, it corresponds to inefficient cooling 
and photo-heating by the UV background. However, note 
that {rih) does not enter the probability distribution (90) 
as long as it is small (the exponential factor in (|6|) can 
be put to unity) which is the case in practice. Thus, our 
results do not depend on the cutoff Tmin- We also get: 



1 



f-^f3, min 



2(3 - P) 



(1 - A^min) 



(93) 



where A*min is the exact minimum value of the magnifica- 
tion obtained in (p2[). As expected, we recover the factor 
(l — fj-min) since the physical process is the same. The term 
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3/^/(2(3 — /?)) of order unity comes from the approxima- 
tions involved in our calculation (assumption of small im- 
pact parameter and approximation (1 + A) ~ A). From 
the results obtained in the previous section, we know that 

/Zniin — 1 so that WC Still havC /i/J^min — 1 and ~ ^ — 1. 

Moreover, the probability distribution ( pO| ) is only valid 
for P{iJ.f}) <C 1 as explained above, that is for ^ 1, 
hence the value of np.jaiu plays no role as long as it re- 
mains small as compared to 1. 

We compare in Figj^ the probability distribution P(/z) 
obtained from ( |90|) for /3 = 1.4 and (3 = 2 with the results 
obtained in Sect!^ from (^). Note that we have defined 
throughout the magnification n by (po|). For large values 
of the convergence k ^ 1 the approximation ( pO[ ) breaks 
down and /i should be obtained from However, since 
in this article we are mostly interested in the regime k 
we always define ^ by (pO|). Thus, /'(/i) can also be under- 
stood as P{k) with the change of variable (pO|). Then, large 
values of the magnification > 3 correspond to k > 1 and 
to strong lensing events with multiple images. We can see 
in Fig.0 that the probability distribution P(/x) from ( pO| ) 
is of the same order as the results from (|4|) in the range 
1.5 < /i < 3. Indeed, we count the same mass and the 
same collapsed halos. Our results are similar to those ob- 
tained by Porciani & Madau (1999) who got a probability 
of a few 10~^ to have a strong lensing event for a source at 
Zs — \ with (3 = 2. However, we can see that (^0|) is quite 
sensitive to the slope (3 of the halos. Indeed, for /? = 2 we 
get P(/i) ^ 11^^ while for (3 = 1.4 we have P{n) ~ ijT^ . On 
the other hand, below the exponential cutoff the prescrip- 
tion developped in Sect.^ leads to P(/i) down to 
P(yLt) ~ . However, since is small this pure power- 
law regime does not really appear as the corrections due to 
the peak at ^ 1 and the exponential cutoff are not neg- 
ligible. Nevertheless, we can clearly see the extended large 
H tail of P(/i) and its non-gaussian behaviour. In partic- 
ular, we recover the trend seen in numerical simulations 
(e.g. Wambsganss et al.l997). As we explained above the 
formulation (^) cannot describe the regime /i ^ 1 while 
for large magnifications it predicts a higher probability 
P(/Lt) since it leads to a power-law tail instead of an expo- 
nential cutoff. This is directly linked to the strong growth 
with p of the p— point correlation functions at small scales 
implied by this model, which does not seem to be com- 
patible with numerical simulations as argued in Valageas 
(1999a). On the other hand, if we only use ( pO| ) to obtain 
the weak lensing effect due to inner galactic halos where 
baryonic matter dominates the density field we would get 
lower probabilities P(/i) since the matter content and size 
of these objects would be smaller, so that up to /i < 3 
the probability distribution (^) would dominate. Thus, 
the formalism developped in Sect.^ is better suited to ob- 
tain the magnification of distant sources by weak lensing 
effects. 

We show in Fig. ^ the mean number of halos {rih) along 
the line of sight up to the redshift Zg of the source, from 




1 1.5 3 3.5 




1 1.5 3 3.5 




Fig. 7. The probability distribution P(/i) of the magni- 
fication of distant sources located at Zg = 0.5 (dashed 
lines) and at = 1 (solid lines). The curves which start 
at /.t = 1 with no falloff at low /i correspond to the for- 
mulation ( |90| ) where the density field is described as a 
collection of virialized halo with a power-law density pro- 
file ( [70| ) with (3 = 1.4 (steeper falloff) oy (3 = 2 (smoother 
falloff). The curves which show a peak at /i < 1 and a 
falloff at lower /i arc from (]45|). 

(H). We use T,„i„ = 3 lO'' K. We can check that {uh) < 1 
which justifies the approximation (pO|). Thus, large mag- 
nifications /i come from the deflection of the light ray by 
a single object. Of course, as explained above this ap- 
proach cannot handle low /i close to /imin- Indeed, these 
lines of sight do not cross totally empty regions but low- 
density patches of matter so that P(/i) for ^ — > /Ltmin, 
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Fig. 8. The mean number of halos {uh) along the line of 
sight up to the redshift Zs of the source. The solid curve is 
for a critical universe, the dashed curve for an open cos- 
mology and the dotted line for a flat low-density universe. 



as in Fig.^, while the formulation 
e~<"''>^D(/^/3) in = 0. 



leads to a Dirac 



6. Derivation of cosmological parameters 

In practice, one uses the Type la supernovae as standard 
candles in order to derive the cosmological parameters Qm 
and Oa from the observed relation redshift <-!■ luminos- 
ity distance (e.g. Perlmutter et al.l999). However, even 
if these sources are perfect candles with no intrinsic dis- 
persion nor instrumental noise, the weak lensing effects 
discussed in the previous sections will introduce some dis- 
persion and some bias as the supernovae will be randomly 
magnified by the density fluctuations located along the 
line of sight. Thus, the luminosity distance measured by 
the observer will differ from the actual one by a small 
deviation which we can relate to the magnification /i by: 



VobsiZs) 



(94) 



where Vobsizs) is the observed distance defined in (|T8|) 
seen as a function of ^m.obs and O.A,obs which are thus 
determined by the observation (if there is no distortion 
by weak lensing: /i = 1, these parameters are equal to 
the actual cosmological parameters VLm and J^a)- Here 
we assume the absolute magnitude of the supernovae is 
known, for instance from local or low redshifts observa- 
tions. We used the fact that the luminosity distance dL{z) 
obeys: dL{z) = (1 + z)'D{z). Thus, for a given redshift Zs 
of the source and a peculiar cosmology (r2,„,i7A) of the 
actual universe we can obtain "confidence regions" in the 
(^m,o6s, ^K,obs) plane for the observed parameters O.m,obs 
and OA.ofjs- This means that for any V with < V < 1, 
the "observed universe" has a probability V to lie within 
the domain TZ-p- Thus, for any point {^m,obsi^K,obs) we 
calculate through (^) the weak lensing magnification 
(jL which is needed so that the observer would measure 



[a 



m^obs 5 



^A.obs) in an actual universe defined by {Qrm ^a)- 



Then, from the intervals of confidence (|6C|), displayed in 
Fig.^, wc obtain the probability V (defined by the con- 
straint that either /i = ^_ or = /i+) that such a magni- 
fication is realized, hence that such a cosmology is derived 
from observations. In this way we obtain confidence re- 
gions TZ-p in the {i^m.obs,^A,obs) plane. Of course, from 
( p4[ ) a measure at a single redshift only provides a value 
for the distance I?ofcs (hence for the deceleration param- 
eter qq at low redshift — s- 0). Hence the parameters 
^m.obs and r^A.obs are only constrained to lie on a line in 
the (Xlm.obsT^A.obs) plane. As a consequence, the regions 
TZ-p are unbounded strips in the {flm,obs,^A,obs) plane. 

We display in Fig.|| these "confidence regions" TZ-p we 
obtain for the three cosmologies we have studied in de- 
tails in the previous sections, for the redshifts Zg = 0.5 
and Zs = 1. At a given redshift Zg, we show the domains 
TZ-p defined hy T' = 68% (dashed region within two solid 
lines) and hy T^ — 95% (within the two dashed lines). This 
means that from one observation at a given redshift the 
observer will conclude with a probability of 68% that the 
cosmological parameters {^m,obs,^A,obs) lie in the dashed 
domain. We also display the boundary line £niin given by 
Mmin within (^) (solid fine on the left side). Thus, be- 
cause of the lower bound /z„iin the observed parameters 
{^m,obs,^A,obs) cannot lie to the left of the line £min (of 
course, here we only consider the effects of weak lensing). 
The cross is the point (fimjf^A)- For a given redshift, all 
these boundary lines are parallel. Moreover, their slope 
(when seen as a function ^A,obs of Om^obs) increases with 
Zg. Indeed, as emphasized by Perlmutter et al.(1997) the 
parameters ftm,obs and flA,obs enter T>obs{z) with different 
powers of (1 -I- z) so that the observed distance T>obs{z) is 
not a function of qobs (except in the limit Zg — > 0) but of 
a combination of flrn,obs and flA,obs which varies with Zg. 
This implies the drift with Zg of the slope of the bound- 
ary lines we defined above. Of course, this is the reason 
why observations can simultaneously constrain Qm and 
flA (provided one has a finite range of source redshifts). 
The sign of the slope of these strips translates the fact 
that for the same flm and redshift a flat universe corre- 
sponds to a larger distance T>{z) than an open geome- 
try. Of course, the strong asymmetry of the probability 
distribution P[^i), and of the intervals leads to 

asymmetric regions TZ-p. Thus, as 7-" increases the strip TZ-p 
grows but it is bounded on the left by £min while it can 
extend to infinity to the right. Indeed, as can be seen from 
( p^ larger /i corresponds to larger ^m,obs and smaller 
T>obs- Thus, the most likely values of (rim.obs, ^A,obs), cor- 
responding to /imax < 1, lie slightly to the left of the point 
(f2„i, ^^a), while for large T' the domain TZ-p shows an ex- 
tended tail towards large ^m.obs- The strips are larger for 
the critical universe which had a higher dispersion for the 
probability distribution P(/x), see Fig.||. We note that al- 
though the spread due to weak lensing cannot make a 
critical universe appear as Vim = 0.3 (due to the cutoff 
/^min) while a low-density universe = 0.3 has a negli- 
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Fig. 9. The "confidence regions" TZ-p for the observed pa- 
rameters ftm,obs and flA,obs, for three actual cosmologies 
and the redshifts Zg = 0.5 and Zg = 1. The dashed region 
within two solid lines corresponds to V = 68%: within 
68% of the cases the cosmology will be observed to lie 
within this domain in the {^m.obs,^A,obs) plane. The re- 
gion within the two dashed lines corresponds toV = 95%. 
The solid line £min on the left corresponds to the lower 
bound /Ltmin: the observed parameters (Jlm,obs, ^A.obs) can- 
not he to the left of this curve. The boundary lines ob- 
tained for a larger redshift are steeper (i.e. closer to the 
vertical) . 

gible probability to appear as = 1, the effect of the 
weak lensing is not negligible. Thus, two observations at 
redshifts z^ = 0.5 and = I only determine f2„i within 



an interval Afi™ ^ 0.3. For instance, for the low-density 
flat universe (lower panel) we have: 



V = 68% : 0.17 < nrn,obs < 0.45 , 0.52 < flA^obs < 0.9 
-P = 95% : < rirn^obs < 0.7 , 0.25 < flA^obs < 1.2 



(95) 



Note that if both observations (at Zg — 0.5 and Zg — 1) 
are independent the intersection of the regions labelled 
V — 68%, used in (p5|), corresponds to a probability 
0.68^ = 46%. These uncertainties are larger for both other 
cosmologies we consider in Fig.^ Of course, by averag- 
ing over many observations at a given redshift one dimin- 
ishes this inaccuracy (one expects that AD,m,obs roughly 
decreases as 1/^/N). Moreover, we can check in the fig- 
ure that observations can unambiguously discriminate be- 
tween fljn — 0.3 and = 1. In the case of a low-density 
universe one can also clearly discriminate between f^A ~ 
and r^A = 1 — ^m- 

In order to compare the effect of weak gravitational 
lensing on the derivation of the cosmological parameters 
with the inaccuracy due to the intrinsic magnitude disper- 
sion <tb of the sources, we show in Fig.RjJ the quantity: 



dM 



B 



'B 



1/2 



(96) 



Here Mb is the absolute magnitude of the supernova and 
^^Biensing the magnitude dispersion due to weak lens- 
ing. Thus, we have: 



InlO dQ.r. 



m,ohs 



2.5 



^=1 



Vinioy 



1/2 



(97) 



where dilm.obs/dfJ, is obtained from (p^. For low-density 
universe we eliminate r^A.ofjs in (|9j) by assuming the right 
cosmology (flat or f^A.obs = 0) so that Vobs is a function 
of ilm,obs alone. In the case of a critical universe we dis- 
play the results we obtain when we assume ^A.obs — or 
^A,obs — ^ — ^m,obs- Wc use for all redshifts ctb — 0.17 mag 
for the lightcurve-width-corrected luminosity dispersion 
of SNela, see Perlmutter et al.(1999). The solid lines in 
Fig.|l^ show Afljn^obs from (|97|). The dashed curves show 
the inaccuracy on ^rn.obs due to aB alone while the dot- 
ted curves show the effect of alone. We can see that 
for low-density universes the error due to the intrinsic dis- 
persion aB dominates up to < 2 while for a critical 
universe the weak lensing contribution already dominates 
for Zg > 1.4. Note also that the inaccuracy of the mea- 
sure of ^m,obs is much larger for a critical universe, which 
has a non-negligible probability to appear as an open uni- 
verse with r2„i — 0.5 or a flat universe with Qrn = 0.6 
for Zg < 1. In particular, note that going to high redshift 
Zg > 1 does not increase the accuracy of the determina- 
tion of by much. Thus, the minimum dispersion of 
^rn.obs is Aflm.obs 0.15 for low-dcnsity universe and 



0.4 for a critical universe. Of course, one can 



reduce the uncertainties by observing many supernovae. 
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Fig. 10. The uncertainty A^jn,obs of the observed param- 
eter rtm from sources at redshift z^, for three cosmologies. 
The solid lines take into account both the intrinsic disper- 
sion as of supernovae and the fluctuation of the mag- 
nification by weak lensing. The dashed curves show the 
influence of ctb alone (case with = 0) while the dotted 
curves represent the influence of weak lensing alone (case 
with as = 0). For an actual universe which is critical 
{Q,n = 1) the lower curves correspond to the assumption 
of a flat universe and the upper curves to f^A = uni- 
verses. For low-density universes we assume the correct 
cosmology (flat or open). 

Note that this analysis does not take into account the 
non-gaussian behaviour of the probability distribution of 
the magnification /i. This was studied in Figjoj Moreover, 
as can be seen in Fig.^ the uncertainty on flm.obs due to 



weak gravitational lensing is larger than the estimate ( |97 
because of the degeneracy in the plane {i^m,obs,^A,obs)- 
Indeed, although observations at two different redshifts 
Zg — 0.5 and Zg = 1 remove this degeneracy the inter- 
section of the domains TZ-p remains elongated along an 
axis roughly parallel to r^A.obs = ^m.obs- Hence the uncer- 
tainty on ilm,obs, given by the projection onto the Qm^obs- 
coordinate of the length of this region along this axis, is 
larger than the value obtained from ( |97| ) shown in FigjlO| 
which corresponds to a cut of the domain TZ-p along the 
axis flA,obs = or QA.obs = 1 - ^m,obs- This is why the 
values obtained in (|95| ) are higher than those one would 
derive from (^7|). 

7. Bias 

The random magnification of distant sources by weak 
lensing induces some bias in any observed sample. For 
instance, the fraction of gravitationally lensed SNela in- 
creases for bright magnitudes due to the falloff of the lumi- 
nosity function of the sources which means that the con- 
tribution from lower luminosity SNela which have been 
amplified by weak lensing becomes more important. We 
assume that the luminosity function of the SNela at red- 
shift z is a gaussian, as a function of the absolute magni- 
tude Mb: 



HMb;z) 



2tt(tb 



(98) 



Here Mbo{z) is the mean absolute magnitude while as 
is the intrinsic magnitude dispersion of the sources. Next, 
we define the bias B{< Mb, z) as the ratio of the number 
of SNela observed with a magnitude brighter than Mb, 
at redshift z, by the number of SNela which actually are 
brighter than this luminosity threshold: 



B{<Mb,z) 



JZ, dM's^M's,z) C(M^-M,,/..= dfiPifi) 



J^'^dM'^mi'B,^) 



where P(/i) is the probability distribution of the magnifi- 
cation. Using (98) we obtain: 



/"OO 

B{< m,z)= / dfiP{fi) 
Jo 



(_ m 



erfc - 



2.5 



log/i 



where we defined the "reduced magnitude" m by: 

Mb - Mbo 
m = 



(99) 



(100) 



and crfc(a::) ~ dt e * is the complementary er- 

ror function. We can see from (691) that: 



(101) 



m — !■ cx) : B{< m,z) ^ 1 
m — !■ — oo : B{< m, z) ^ oo 
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n„-o.3 , n,-o 



Fig. 11. The bias B{< m, 
and Zs = 1 (solid lines). 



Indeed, for m ^ oo one recovers all SNela (the survey 
is complete) while for m ^ — cxd the local slope of the 
luminosity function is increasingly steep so that the effect 
of weak lensing gets larger. 

We present in Fig.^our results for the source redshifts 
Zs = 0.5 and Zg = 1. We can see that the bias is larger 
at higher redshift. Indeed, as can be seen from Fig.^ the 
large /x tail of the probability distribution P{^J^) is more 
important at higher z because of the larger variance 
This implies that the effect of weak lensing is higher (it is 
more likely for SNela to be "amplified" by density fluctua- 
tions) hence B{< m, Zg) is larger. Note that the bias is not 
negligible: we already get B{< m,Zs) ~ 2 for m = — 2 at 
Zs = 1, which only corresponds to a "2 sigma" deviation 



from the mean Mbo for the magnitude of the observed 
supernova. 

These distortions due to weak lensing lead to an ap- 
parent luminosity function $l of the sources: 

^l{Mb;z)^ dfiP{fi)^MB + 2.5\ogfi;z) (102) 
Jo 

Thus, the mean magnitude (Mb) Msth of ^ survey limited 
by the absolute magnitude threshold AIbih is: 



B/AlBth 



M 



1 



BO 



/d/i^(M)crfc(-I^H^^^±^^) 




(Tb exp 



^thCTB + 2.5 log /i 



V2< 



CfB 



-2.5 log/i erfc — 



mthCTB + 2.5 log 



O-B 



(103) 

The first term (with the exponential) mainly corresponds 
to the fact that the survey is limited by the upper mag- 
nitude MBth which implies that {MB)MBth < Mbui- It 
vanishes for Mbih — > oo when the survey is complete. The 
second term (with the prefactor 2.5 log /i) does not go to 
zero even if Mbih oo: it is due to the distortion of the 
luminosity function which implies that the mean of is 
no longer Mbo- Indeed, we obtain at the lowest order in 



{Mb)oo = Mbo -2.5 J dfi P{ii) log/x 



2.5 . 

^^^''+2ha0^^ 



(104) 



5) for Zs ~ 0.5 (dashed lines) and' 




(105) 



V\/21nl0 



As we can check in (105) the dispersion due to weak lens- 
ing increases the bright and faint tails of the luminosity 
function. However, this distortion is not symmetric, even 
at the lowest order in because (log^) < 0. Hence the 
faint part of ^{Mb) shows a larger increase than the very 
bright part and {Mb)qo > Mbo- As can be seen from (105) 
and Figjl^ the effect of weak lensing is more important if 
the survey is not complete, hence the dependence on 
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n„ = o.3 . n. = o 



n„ = o.3 . n. = o.v 



Fig. 12. The mean magnitude {MB)MBth of ^ survey 
limited by the apparent magnitude Msth for = 0.5 
(dashed Hnes) and Zg = 1 (solid hues). We note nith = 
{Msth — MBo)/crB the deviation of the threshold Msth 
from the intrinsic mean Mbo in units of the dispersion 

(7B- 

(i.e. on the redshift Zg of the source) is larger. In particu- 
lar, we have: 

2 

rr D 

mth < -1 , < 



2.5m,th 



(106) 



We present in Fig.|l^ the deviation from the intrinsic 
mean Mbq of the average magnitude measured by a sur- 



vey with the upper absolute magnitude threshold Msth, 
given by the "reduced magnitude" threshold mth as de- 
fined in ( |100| ). For bright magnitude thresholds the ap- 
parent mean is close to Msth while for faint Mbih the 
survey is almost complete and {MB)MBth — ^^bo- As ex- 
plained above, even for MBth — > oo the observed mean is 



not equal to Mbo- However, as we can see in (104) and 
in Figjl^ the deviation is very small for the redshifts of 
interest Zs ^ 1 since the variance of the magnification 
is small, see Fig.^ In particular, we note that the curve 
{MB){MBth) shows a very small dependence on redshift. 
In practice, in order to measure the cosmological param- 
eters and l^A one observes SNela at low {zg ~ 0.04) 
and high {zg ~ 0.8) redshift. The low Zg data gives the 
normalization of the curve mBapp{z), where niBapp is the 
apparent magnitude, while the large Zg data constrains 
the cosmological parameters. Thus, the determination of 
the cosmological parameters is only sensitive to the differ- 
ence between the deviations (Mb) MBth ~ ^bo at low and 
large Zg. Provided both surveys have a sufficiently faint 
absolute magnitude threshold: MBth > Mbo + 2crB, the 
bias AAIb due to weak lensing effects is: 



AMb{zi,Z2) = {{MB)MBtu2 - Mbo){z2) 

-{{Mb)mbu.,-Mbo){zi) 



(107) 



where zi < Z2 are the redshifts of both surveys. Thus, if 
both surveys are almost complete we have: 



AAfB(zi,Z2) < 



2.5 

2 InlO 



(108) 
Fig.0 



For 22 < 1 we get AMb(zi,z2) < 5.4 10"^, see 
Thus this effect is negligible for the determination of cos- 
mological parameters by SNela. Of course, if the low and 
large redshift surveys have bright magnitude thresholds 
MBth ^ Mbo which are different it is not possible to es- 
timate rtm nor r^A- This does not seem to be the case in 
practice (Perlmutter et al.l999). 

In order to derive from observations the cosmologi- 
cal parameters (f2m,^^A) one draws a redshift «-> appar- 
ent magnitude diagram. Moreover, the apparent magni- 
tude ruapp is "corrected" thanks to the light-curve width- 
luminosity relation (e.g. Perlmutter ct al.l999). Of course, 
because there is some scatter in this latter relation there 
is still a small lightcurve- width-corrected dispersion aB — 
0.17 mag. The curve mapp{z) depends on the cosmology, 
through the luminosity distance di, which allows one to 
derive (17™, JIa): 



Mp, 



(109) 



where the observed apparent magnitude has been cor- 
rected for K and extinction corrections. However, due 
to the apparent magnitude threshold of SNela surveys a 
Malmquist bias appears at large redshifts where the mean 



20 



P. Valageas: Weak gravitational lensing effects on tlie determination of f2,„ and Qa from SNela 



O.S 0.4 0.6 O.B 



n„=0-3 , n, = o 



0.2 0.4 O.S O.B 



n„=0-3 , v.. = 0-v 



Fig. 13. The redshift ^ apparent magnitude diagram for 
three cosmologies. We display the quantity {mapp)—MBo+ 
5 log(/i) as a function of the redshift Zg of the source. The 
dotted curves show the apparent magnitude-redshift re- 



lation (109) obtained for a survey which is not limited 
by a flux threshold. The solid (resp. dashed) curves show 
the effect of an apparent magnitude threshold when weak 
gravitational lensing is (resp. is not) taken into account. 

apparent magnitude of supernovae is close to this thresh- 
old. Indeed, the curve one obtains from observations is not 
(|lO|) but: 

dL 



{mapp)MstH(z) = {MB)MBtH(z) + 5 log 



lOpc 



(110) 



where (Mb) Msthiz) mean magnitude of a survey 

limited by the absolute magnitude threshold MBth{z), as 



in (103). The latter is obtained from a given apparent 
magnitude threshold mappth by: 



mappth = MBth{z) + Slog 



lOpc 



(111) 



We display in Figjl^ the curves {niapp) {zs)—Mbo+5 log(/i) 
we obtain for three cosmologies (the term \og{h) removes 
the dependence on the Hubble constant). The dotted 
curves, which correspond to ( |109| ) (no magnitude thresh- 
old), show the (small) dependence on cosmology of the ap- 
parent magnitude-redshift relation: low-density universes 
lead to larger luminosity distances hence to larger appar- 
ent magnitude (fainter object) at fixed Zg, this effect is 
larger for the flat model than for the open case. The solid 
(resp. dashed) curve shows the effect of an apparent mag- 
nitude threshold when weak gravitational lensing is (resp. 
is not) taken into account. The threshold is chosen so that 
MBth{z) — Mbq at redshift zth — 0.7. Thus, we see in the 
figure that at lower redshifts Zg < 0.6 the apparent mag- 
nitude threshold of the survey plays no role: all curves su- 
perpose onto ( 109 ). At higher redshift Zg > 0.6, because of 
the threshold ruappth, the survey only detects the bright- 
est SNela which means that the average {"mapp) Msthizs) 
is biased towards small magnitudes (large luminosities). 
This leads to a clear de viat ion of the observed magnitude- 
redshift relation from ( |109| ). Thus, the break in the curve 
clearly marks the redshift beyond which the cosmologi- 
cal parameters cannot be derived from observations with 
this apparent magnitude threshold. Taking into account 
the weak lensing effects (solid lines) slightly amplifies this 
bias towards large luminosities because the random mag- 
nification by density fluctuations along the line of sight 
increases the large apparent luminosity tail of the SNela 
distribution. 

8. Conclusion 

In this article, we have shown how one can obtain the prob- 
ability distribution of the magnification of distant sources 
by weak gravitational lensing, using a realistic description 
of the density field which has already been checked against 
numerical simulations of structure formation within hier- 
archical scenarios. Thus, this work improves the results 
obtained by previous analytical studies which used for in- 
stance a "Swiss cheese" model to describe the universe 
(Kantowski 1998) or considered a collection of virialized 
halos amid an empty space. We recover the behaviour ob- 
served in numerical simulations, which is not surprising 
since we consider similar density fields. Thus, the proba- 
bility distribution of the magnification shows a maximum 
at a value slightly smaller than the mean (/i) = 1 and 
it shows an extended large fj, tail. Moreover, the vari- 
ance increases at larger redshifts while the deviation 
from a gaussian gets higher at lower redshifts. The advan- 
tage of our approach is that we obtain a direct connection 
of the weak lensing properties with the characteristics of 
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the underlying non-linear density field. In particular, the 
non-gaussian behaviour of the magnification is expressed 
in terms of the non-gaussian properties of the density field, 
through its many-body correlation functions. 

Then, we have applied our results to the magnification 
of distant Type la supernovae. We have shown that the 
inaccuracy introduced by weak lensing in the derivation of 
the cosmological parameters is not negligible: AVlm ^ 0.3 
for two observations at Zs — 0.5 and Zs = 1. However, ob- 
servations can unambiguously discriminate between flm = 
0.3 and ilm = 1- Moreover, in the case of a low-density 
universe one can clearly discriminate between JIa = and 
fl\ = 1 — rirn- Besides, the accuracy increases as the num- 
ber of SNela gets larger (there are already 42 available 
SNela, see Perlmutter et al.l999). On the other hand, if 
it were possible to measure the distortions due to weak 
lensing one would obtain some valuable information on 
the properties of the underlying non-Hnear density field, 
since we have shown that the probability distribution of 
the magnification can be directly expressed in terms of 
the probability distribution of the density contrast at the 
non-linear scale (typical of present galaxies) where the lo- 
cal slope of the initial linear power-spectrum is n = —2. 
However, this would require a rather high accuracy of the 
observations, as we have shown that the probability dis- 
tribution of the magnification is not very extended (the 
typical deviation is of order A/i ~ 0.08 at Zg = 1, hence 
Amag ~ 0.08). A more detailed discussion of the proper- 
ties of the p.d.f. P{k) and P(Map) of the convergence k 
and the aperture mass M^p is presented in other articles 
(Valageas 1999b; Valageas 2000; Bernardeau & Valageas 
2000). 
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